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The orbitally degenerate A-site spinel compound FeSc2S4 has been experimentally identified as a "spin- 
orbital liquid", with strong fluctuations of both spins and orbitals. Assuming that the second neighbor spin 
exchange J2 is the dominant one, we argued in a recent theoretical study [Chen et al, Phys. Rev. Lett. 102, 
096406 (2009)] that FeSc2S4 is in a local "spin-orbital singlet" state driven by spin orbit coupling, close to a 
quantum critical point, which separates the "spin-orbital singlet" phase from a magnetically and orbitally or- 
dered phase. In this paper, we refine further and develop this theory of FeSc2S4. First, we show that inclusion of 
a small first neighbor exchange Jj induces a narrow region of incommensurate phase near the quantum critical 
point. Next, we derive the phase diagram in the presence of an external magnetic field B, and show that the 
latter suppresses the ordered phase. Lastly, we compute the field dependent dynamical magnetic susceptibil- 
ity x(k, w; B), from which we extract a variety of physical quantities. Comparison with and suggestions for 
experiment are discussed. 

PACS numbers: 71.70.Ej,71.70.Gm,75.10.-b,75.40.-s 



I. INTRODUCTION 

Among all the magnetic spinelsl^™ ^ ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 ^^ A- 
site spinel FeSc2S4 is particularly intriguing in that its frus- 
tration paramete J 14 ! 15 ! / > 1000 is one of the largest ever 
reported Indeed, even though the material clearly exhibits 
well-formed local moments interacting with a characteristic 
energy given by the Curie-Weiss temperature |0cw| = 45K, 
no sign of magnetic ordering has been found down to the low- 
est measurable temperature of 50mK. Moreover, FeSc2S 4 is 
interesting because it has not only spin but also orbital degen- 
eracy. The Fe 2+ ion at the A-site is in a 3d 6 configuration, 
whose five-fold degeneracy is split by the tetrahedral crystal 
field into a lower e g doublet and an upper ti g triplet. The six 
electrons in the 3d shell are Hund's rule coupled, yielding a 
high spin configuration with S = 2 and a two-fold orbital de- 
generacy due to a hole in the lower e g doublet. Besides the 
five-fold spin degeneracy and the lattice degrees of freedom, 
the orbital degeneracy gives an additional contribution to the 
specific heat. This has been confirmed experimentally.™ 

Commonly, orbital degeneracies are relieved by a Jahn- 
Teller type structural distortion, that leads to orbital order at 
low temperatures. However, in FeSc2S4 no such distortion has 
been observed. Hence, both the spins and orbitals remain frus- 
trated and continue to fluctuate down to the lowest measured 
temperature, a situation for which the term "spin orbital liq- 
uid" (S OL) w as coined. A SOL state was also suggested for 
L1N1O2 ESIIZI But this proposal has been questioned recently 
in Ref . 1 8 , where it was argued that the unusual behavior of 
LiNi02 is due to disorder effects. Therefore FeSc2S4 remains 
as the best candidate for a SOL. 

In a undistorted lattice, exchange interactions or spin-orbit 
coupling can split the orbital and spin degeneracies of a sin- 
gle Fe 2+ ion. To study these possibilities we introduced in 
Ref. [19] a model that contains a "Kugel-Khomskii"-type 2() 
spin-orbital exchange interaction as well as the lowest-order 
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FIG. 1: (color online). Schematic phase diagram as a function of 
temperature T and x, the ratio of exchange to spin-orbit interac- 
tion. A quantum critical point (QCP) at x c separates the spin-orbital- 
singlet (SOS) state from the ordered phase. Within the order phase, 
at x C 2, there is a commensurate-incommensurate phase transition. 
"QC" denotes the quantum critical regime, "CO" the commensu- 
rate antiferromagnet with orbital order, and "IC" the incommensurate 
spin and orbital order. Inclusion of a small first neighbor exchange 
Ji induces a second region of incommensurate phase (shaded area). 



symmetry-allowed atomic spin-orbit coupling. The exchange 
interactions favor spin and orbital order, whereas the on-site 
spin-orbit coupling leads to the formation of a local "spin- 
orbital singlet" (SOS). In FeSc2S4 there is strong competition 
between these two interactions. We argued in Ref. [19] that 
FeSc2S4 is in the SOS state, close to a quantum critical point 
(QCP) that separates the SOS phase from the ordered phase 
(see Fig. [TJ. This QCP seems to be rather analogous to that 
appearing in spin-dimer materiald^l or bilayer Heisenberg 
models, 23 with the two orbital states playing the roles of the 
two members of a spin dimer, or the two bilayer states. How- 
ever, we will discover at least one surprising difference below. 

Previously, we have argued that, to a first approximation, 
FeSc2S4 can be described by a simplified " J2 — A modei"p2l 
which only contains on-site spin-orbit and next-nearest neigh- 
bor (NNN) spin exchange interactions. Within this approxi- 
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FIG. 2: (color online). Phase diagram of the J2 — A model, Eq. {14) , 
as a function of magnetic field B and coupling ration x. In zero 
magnetic field the phase transition between the spin-orbital singlet 
(SOS) phase and the commensurate ordered (CO) state occurs at 
x c = 1/16. 



mation the two fee sublattices of the A-site diamond lattice 
decouple completely. The aim of the present paper is to refine 
the model of Ref. [19] and to calculate more detailed physical 
properties for comparison with experiments. 

We summarize the results here. We first consider the ef- 
fects of weak intersublattice nearest-neighbor (NN) exchange 
J\. We find that it induces a narrow range of incommensu- 
rate magnetically ordered phase in the vicinity of the QCP, on 
the ordered side (see Fig. [TJ. Next, we consider the effects 
of an external magnetic field B on the QCP in the J 2 — A 
model. Quite surprisingly, we find that the magnetic field de- 
stroys the spin order and the QCP shifts from x c — 1/16 to 
a larger value (see Fig. EJ. This is exactly the opposite trend 
to that observed in the spin-dimer and bilayer models. Within 
mean field theory we calculate both the uniform and staggered 
magnetization as a function of B. Finally, we use the random 
phase approximation (RPA) to compute the field-dependent 
dynamical spin susceptibility Imx(fe, ui, B) in the SOS phase. 
From an analysis of the pole structure of x(fe,o;) we obtain 
the dispersion of the low-energy collective modes. Consistent 
with the effect of the magnetic field on the phase diagram, we 
find that the gap in the SOS phase increases with B. 

The remainder of the paper is organized as follows. In 
Sec. [TTJ we define the model Hamiltonian and summarize the 
principle results of Ref. [19] The weak inter-sublattice inter- 
action is introduced in Sec. [In] where we derive the result- 
ing changes to the phase diagram. 
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we return to 



In Sec. 

the minimal NNN model, but include the effects of an exter- 
nal magnetic field. The inelastic structure factor (dynamical 
spin susceptibility) is calculated by the RPA approximation 
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Sec 
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V] We conclude with a summary and discussion in 
Some technicalities are given in the appendix. 



H. MODEL DEFINITION 

The Hamiltonian contains two terms: an exchange interac- 
tion H ex and an atomic spin-orbit coupling H l , 



(la) 



The term 7i ex describes spin and orbital exchange interactions 
as well as couplings between spin and orbital degrees of free- 
dom. Using microscopic considerations and symmetry con- 
straints one can show that 7i ex is of the forrrU^ 



[Jij S, • S, + 2K ij T i • T,( 1 + S, • S,); .(lb) 



where Si is the S — 2 spin operator at site i. The orbital 
degrees of freedom are described by T = 1/2 pseudospin op- 
erators T, that act on the x 2 — y 2 and 3z 2 — r 2 orbitals in the 
e g subspace. An analysis of the exchange paths linking two 
A-sites shows that both first and second neighbor exchange 
pa ths a re of comparable length and have similar multiplic- 
itytZlSU rp n j s SU gg eS |; S that there is a substantial NNN interac- 
tion and it is therefore necessary to keep in Eq. ( [TbJ the sum 
over both the NN and NNN sites. For convenience we set 
•hi = Ji or J2 when ij are first and second neighbor sites, 
respectively (and similarly for 

The second term in Eq. ( lai, 7Yq, is the on-site spin-orbit 
coupling which arises from second order perturbation theory 25 



n l 



{V3T?[(Sf f - (Sf) 2 } + Tf[3(Sf) 2 - S 2 ]} , 



(lc) 



where the coefficient A is estimated from microscopic atomic 
calculations to be A ~ 6Ag/A te . Here, Ao denotes the atomic 
spin-orbit interaction and A te is the energy separation between 
the e g and t^g states. Note, that 7Yq immediately leads to a 
splitting of the ionic degeneracies because the ground state of 
Hq is the non-degenerate spin-orbital singlet 



1 

V2 



x 2 - y 2 ) \S Z = 0) + 



) IS 2 



(2) 



i|3z 2 



\S Z 



-2) 



Whether FeSc2S4 is in such a SOS phase depends on the ratio 
between exchange and on-site spin-orbit interactions 



max { Ji, J 2 , Ki, K2} /A. 



(3) 



We note that we have not included spin-orbit effects such as 
Dzyaloshinskii-Moriya interactions in the exchange Hamilto- 



nian, Eq. ( lb 1. This may be surprising since the on-site spin- 
orbit coupling in Eq. ( |Tc] > plays a crucial role in our analysis. 
However, we expect that the spin-orbit corrections to the ex- 
change are smaller than the leading order exchange couplings 
by a factor of order Ao/A te . This makes them subdominant 
both to the isotropic exchange couplings in Eq. ( |Tb] > and the 
on-site spin-orbit interaction in Eq. ( fTc] i, which are compara- 
ble in FeSc2S4 (this indeed defines the location of the QCP - 
see below). 

From a comparison with experiments we showed in Ref. [19] 
that Ji is antiferromagnetic and the largest among all the 
exchange coupling constants, whereas A is competitive, but 
slightly larger than the NNN spin exchange interaction J2. 
This observation leads us to consider a "minimal" version of 
model ([TJ, which only includes the NNN spin exchange J 2 
and the onsite spin-orbit interaction A (i.e., J\ = K\ = K2 — 
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0). Within this J 2 — A model we demonstrated that FeSc2S4 respectively, 
is in the SOS phase close to the QCP of FigJTJ 

The full phase diagram of Hamiltonian (IT 1 as a function of 
temperature and ratio x, Eq. ([3]), includes a commensurate- 
incommensurate transition within the ordered phase (see 
Fig. [TJ). Deep in the ordered phase when the exchange inter- 
actions are dominant (x 3> 1) the incommensurate spin and 
orbital order (IC) is generally favored by the exchange Hamil- 
tonian ( Taj. With the inclusion of weak spin-orbit interaction, 
Eq. ([lb 1, i.e., with decreasing x, the spin and orbital order be- 
comes commensurate with the spins and orbitals both forming 
a spiral with wavevectors p = 2ir(l, 0, 0) and q = (0, 0, 0), 



III. EFFECTS OF WEAK INTER-SUBLATTICE SPIN 
EXCHANGE INTERACTION 



As mentioned in the introduction, the inclusion of a small 
NN interaction Ji induces a narrow region of incommensurate 
order near the QCP (shaded area in Fig. [TJ. In this section we 
give a derivation of this results using a Landau expansion of 
the effective action. 



First, we note that the Hamiltonian in Eq. ( fT4| i has independent cubic "internal" spin symmetry and cubic "external" space 
group symmetry. We therefore have the symmetry allowed free energy for two decoupled fee sublattices 19 near the QCP, 



MM) = E [^l^,a| 2 + «2El^a| 2 +r|^, Q | 2 ] + 5l ^(|^ a 



ft , a 



fj,=A,B;a b^a 

J2 Sym[ 53 (^, 1 ) 2 (^. 2 ) 2 + <? 4 (^,i) 2 (^, 2 ) 2 + a^lMA^l, 



(4) 



where the order parameters rp^ a are the (real) staggered magnetizations introduced by 



(S,) = 



J2 ~<PaA~1) 2x * « e A sublattice 

a— x,y ,z 

V'B,a(-l) 2;E *' i € B sublattice 

a— x,y ,z 



(5) 



and xf are the canonical half-integer coordinates of the fee lattice with the cubic supercell having unit length. For our con- 
venience we choose the same set of fee coordinates for the two sublattices. In Eq. Q, "Sym" indicates symmetrization with 
respect to both wavevector (lower) and spin (upper) indices. 
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We now introduce weak inter-sublattice interactions. As 
a result, extra terms which couple two fee sublattices will ap- 
pear in Eq. Q. By the symmetry analysis, we obtain the inter- 
sublattice terms 



(6) 



where we have kept the leading quadratic term (others have 
more derivatives) and the most important quartic term, which 
is when inter-sublattice couplings first appear without deriva- 
tives. A number of other quartic couplings involving the two 
sublattices are also allowed, but do not play any important role 
in what follows. 

We now proceed to analyze the Landau theory of Eqs. Q 
and (|6j. First consider the behavior on approaching the QCP 
from the SOS phase. The first instability of the SOS phase 
is signalled by the quadratic part of the action alone, i.e. the 
vanishing of the lowest eigenvalue of the associated quadratic 
form. Due to the linear derivative term, the unstable eigenvec- 
tors are non-constant fields of the form 



tjj b a cos(Sx a 
ip b a sin(5 x a ■ 



i,b), 



(7) 



with S = 7/(4?^), and V'ai^a-fc arbitrary constants. There is 
one linearly independent unstable eigenvector for each a and 
b. The particular form of superposition which is favored in the 
ordered state is determined by the quartic terms. We expect on 
physical grounds that the ordered states will be of spiral type, 
with approximately constant magnitude of spin expectation 
values, and with a "single q" structure. The latter condition 
means that ipA,a is non-zero only for a single a = x,y or z. 
A sufficient condition for a single q structure to be favored 
is that 33,54,55 > (though this condition can be relaxed 
somewhat). In this case, we have 



1pA,a = ^oRe 



(ei + ie 2 )e 



i(5x a +6) 



(8) 



tpB,a = V'oRc (-e 2 + iei)e 



i(5x a +9) 



This corresponds to an incommensurate spiral state of spins 
with the wavevector 



(2tt±5)(1,0 ! 0) . 



(9) 



The unit vectors ei/ 2 define the plane in which the spins ro- 
tate. At the quadratic level, this plane is arbitrary, but it will 
be selected by the quartic terms in the free energy. Presuming 
that the system has axial cubic anisotropy (preferring spins 
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aligned with the x,y or z axes), a (100) or symmetry -related 
plane will be chosen. This is controlled in the free energy by 
the coefficient g 2 , which should be negative to mimic axial 
cubic anisotropy. 

Now consider the evolution of the spin configuration as the 
system becomes more strongly ordered. As \rp \ increases, we 
expect the quartic terms in the free energy to become more im- 
portant, favoring commensurate states in which the spins are 
aligned with the principle axes. To analyze how this occurs, 
we presume that the spins remain in a "single q" structure (a 
spiral rather than a more exotic "spin lattice"), so that i/v,a 
is non-zero only for one a. Furthermore, we assume that the 
spins remain in a single (100) plane. Up to symmetry -related 
choices we take a = x = 1, and ei = x, e 2 = y in Eq. ^Jj. 
Finally, we presume that the fields depend only upon x, as can 
be easily verified is true for the minimum free energy config- 
urations. 

Inserting this into Eqs. Q and d6l, and assuming ipo is con- 
stant, we obtain the reduced free energy density, neglecting an 
additive constant 

/ = ^(d x 0) 2 + a cos 4(0+ fo), (10) 

with k — AvfipQ, a = <72"0o/2. This is a standard sine- 
Gordon model, with incommensuration S. It describes a com- 
petition between a commensurate state in which 9 + 5x is 
constant, and an incommensurate one in which d x 9 < 5 on 
average. The transition between the two states is known as a 
Commensurate-Incommensurate-Transition (CIT), and its lo- 
cation is determined by the condition that the energy of a sin- 
gle domain wall excitation of the commensurate state (a "soli- 
ton") vanishes. In this way we can precisely determine the 
location of the CIT. We find that the CIT occurs foPl 
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FIG. 3: (color online), (a) Uniform magnetization m versus mag- 
netic field B (in units of the spin exchange J2) for different coupling 
ratios x (from top to bottom x = 0.15,0.10,0.05,0.02, respec- 
tively). The arrows indicate the position of the QCP. (b) Staggered 
magnetization n as a function of B (in units of the spin exchange J2) 
for x = 0.15, 0.10, 0.07, 0.063 (from top to bottom curve). 



for J 1 <C J 2 . As expected, for small J\, the incommensurate 
phase studied here is narrow and located only near the QCP 
(see Fig. [TJ. As discussed in Ref. [19] a different incommen- 
surate phase arises for much larger J 2 , when the spin-orbit 
interaction plays a minimal role. This transition to the sec- 
ond incommensurate phase occurs at x c2 ~ 0.61( J2/J1) 2 , 
far from the QCP. 



When 5 < 5 C , the system is in a commensurate state. 

We must now translate this condition back to obtain the 
critical point in terms of microscopic parameters x, x c , and 
J1/J2. In Appendix |A[ we determine the necessary coef- 
ficients in Eqs. Q and |6]) (v 2 , r, <?i, <?2, 7) by deriving the 
effective action by standard Hubbard-Stratonovich methods 
from the microscopic spin-orbital Hamiltonian. The results 
are given in Eqs. ( |A8| > and ( |A9] >. In addition, we require the 
amplitude ipo to minimize the free energy. With this we can 
calculate k and a, and hence solve Eq. (JTTJ for the location 
of the CIT. To get the amplitude, we recognize that at the CIT 
the solution is commensurate, i.e. has the form of Eq. ([8]) with 
6 = 9 = 0. Therefore, we may simply evaluate the free en- 
ergy in Eq. |4]) using this form of the order parameter, and 
minimize over tpo. The result is that 



2(5i+52)' 



(12) 



where at the CIT of course r < 0. Combining the above 
results, we obtain the CIT transition point 
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IV. J 2 - A MODEL IN A MAGNETIC FIELD 

In this section we study the minimal J 2 — A model in an 
external magnetic field B both in the CO and the SOS phase 
(see Fig. [TJ. For definitiveness we assume that the spins in 
the CO phase align themselves along the x axis, whereas the 
magnetic field is applied along the z axis. We use mean field 
theory to calculate the uniform and staggered magnetizations 
and derive therefrom the magnetic phase diagram of the J 2 — \ 
model. 

The minimal J 2 — A model in an external magnetic field 
B contains only onsite and NNN interactions. Therefore the 
diamond lattice decomposes into two fee sublattices with J 2 
playing the role of the NN exchange interaction within each 
sublattice. The Hamiltonian is 



(ij) 



(14) 



where (ij) represents nearest neighbor sites on an fee sub- 
lattice and Hq is defined by Eq. ( [Tcj l. In the absence of a 
magnetic field and when J2/A > x c = 1/16 the system is in 
the CO phase 19 with the spins aligned antiferromagnetically 
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along one of the three cubic axes (here taken to be the x axis). 
To decouple the exchange interaction in Eq. ( |T4] > we employ 
mean field theory with the following ansatz for the average 
spin at site i 



ncos (p • Ti)x 



mz 



(15) 



Here, m and n denote uniform and staggered magnetizations, 
respectively, and are the usual half-integer coordinates of 
the fee sites. In the CO state the spiral momentum takes the 
form p = 2%(1, 0, 0) thereby encoding the antiferromagnetic 
order along the x axis. In the disordered SOS phase the stag- 
gered moment is vanishing, n = 0. With this, the resulting 
single-site mean field Hamiltonian reads 



7if F = h z Sf + h x Sf + W , 



(16a) 



where we have introduced the two effective magnetic fields 



12 J 2 m + B, h x 



-4J 2 n 



(16b) 



At zero temperature the self-consistent mean-field equations 
for the uniform and staggered magnetizations are given by 
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de(h x ,h z ) 

dh z 
de(h x ,h z ) 



(17a) 
(17b) 



with 



x c (B) 



1 

16 



B 2 



16384 Jf 



(21) 



where we retained only the lowest order term in B. The ap- 
proximate result for the uniform magnetization, Eq. ( |19) , de- 
scribes the linear dependence on magnetic field in the regime 
where B is small compared to J 2 (cf. Fig. [3]). It is interest- 
ing to note that this behavior agrees with the measured low- 
temperature magnetic susceptibility in FeSc2S4. Indeed, it is 
found in Ref. 7 that the magnetic susceptibility in FeSc2S4 at 
T — > saturates to a constant value, independent of B. 

To conclude, we find that an external magnetic field leads 
to a suppression of spin ordering. On the ordered side of the 
QCP there is a phase transition to the disordered SOS phase 
as the magnetic field is increased. This behavior is quite dif- 
ferent from a spin singlet phase in a typical spin-only system 
such as TICUCI3 21 or BaCuSi 2 Og 22 , where only the magnetic 
triplet excited states - magnons - respond to the magnetic 
field, whereas the non-magnetic singlet ground state is unaf- 
fected. Hence, the field stabilizes spin order by Bose-Einstein 
condensation of magnons. Here, however, the strong spin- 
orbit interaction leads to very different physics. Specifically, 
SOS is not a spin singlet, but a highly entangled quantum state 
of spin and orbital degrees of freedom. As a consequence, it 
responds strongly to the applied field, and indeed takes better 
advantage of the field than does the ordered Neel state. 



where e(h x , h z ) denotes the ground state energy of the mean- 
field Hamiltonian ([16]). The numerical solutions to these equa- 
tions are presented in Fig. [3] We find that with increasing field 
the staggered magnetization is suppressed, and eventually the 
magnetic order is destroyed. Hence, the critical coupling ratio 
x c moves from x c = 1/16 to larger values with increasing 
field (see Fig. [2j. The uniform magnetization shows a small 
"shoulder" at the critical magnetic field when the staggered 
magnetization vanishes. 

In the neighborhood of the QCP (\x/x c — 1| -C 1) and for 
B small compared to J 2 it is legitimate to expand the ground 
state energy e(h x , h z ) in the effective magnetic fields h x and 
h z . Up to fourth order we have 



e(h x , h z ) 



-2A - - (h 2 x 



2 

A^ 



Ahth 



2.2 



hi). 
(18) 



Neglecting terms of order h z and higher, we obtain the fol- 
lowing expressions for the uniform magnetization 
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SOS, 
CO, 



(19) 



in the SOS and CO phase, respectively. Similarly, the stag- 
gered magnetization is given by 



0, SOS, 

8y/2(x-x c (B)), CO, 



(20) 



V. DYNAMICAL SPIN SUSCEPTIBILITY AND ENERGY 
GAPS 

In this section, we compute the dynamical spin susceptibil- 
ity in the SOS phase with the exchange coupling J 2 treated 
within the RPA. This could be compared to inelastic neutron 
scattering data. The energy gaps to the low-lying collective 
modes are derived from an analysis of the pole structure of 
the dynamical spin susceptibility. 

In the SOS phase we can approximate the full magnetic 
susceptibility of the J 2 — A model x MI/ (k, w; B) by the RPA 
in terms of the non-interaction susceptibility X n U (k , u>, B) of 
the on-site Hamiltonian Hq + BSf [see Eq. (Id]. For the 
xz-component we have 



X XX (K^;B) 



xiTM 
i-j(k) x g x H 



(22) 



where J{k) is the Fourier transform of the exchange cou- 
pling, 



J-(k) = Yj ^cos(k- A), 

{A} 



(23) 



with {A} denoting the 12 NNN lattice vectors. The single- 
site spin susceptibility ~x^ x (u>)B) can be constructed from 
the spectral representation. At zero temperature Xo'i.^B) 
is given by 



XX 

AO 



\(o\s x \j)\ 2 , |<o|s*|j>| s 



iT 



iT 



(24) 
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FIG. 4: (color online) Imaginary part of the angular-averaged 
RPA spin susceptibility in zero magnetic field Im x^y (k, 0) with 
J 2 /A = 0.05 and a damping F/\ = 0.2. 



FIG. 5: (color online). Energy gap A x (in units of the spin-orbit 
coupling A) as a function of magnetic field (in units of the spin- 
exchange J2) for different coupling rations x (from top to bottom 
x = 0.05, 0.06, 0.062, 1/16). 



where |0) and \j) are the ground state and excited states of 
the on-site term H.q + BSf, respectively. The energy dif- 
ference between the excitated state |j) and the ground state 
1 0) is denoted by €j. Finite lifetime effects are parametrized 
by a phenomenological damping T > 0. In order to fa- 
cilitate a direct comparison with neutron scattering data on 
polycrystalline samples we perform a numerical average of 
X^k, ui; B) over the angular components of the wavevector 
k. For a given wavevector magnitude k = \k\ we define the 
angular-averaged spin susceptibility by 



sin ( 



X xx (k,co;B), 



(25) 



where 8 and (j> describe the direction of the wavevector 
fe. Since the inelastic neutron scattering intensity is pro- 
portional to the imaginary part of xtvi^i B) we com- 
pute lmxfy(k, ui; B) as a function of energy transfer ui and 
wavevector magnitude k. Fig.|4]displays the numerically cal- 
culated dynamical spin susceptibility lmx xx (k, ui; 0) in zero 
magnetic field. The excitation minima near k = 2ir/a and 
k = 5ir/a agree well with the neutron scattering data on poly- 
crystalline samplesP 

The dispersing excitation branch shown in Fig.[4]is a collec- 
tive mode associated with zeros of the real part of the denom- 
inator in Eq. ( |22) . For the unaveraged dynamical spin sus- 
ceptibility lmx xx (k, u>, B) the frequency minimum A x of the 
dispersing collective mode occurs at k = 27r(l,0,0). Upon 
approaching the QCP from the disordered side and for small 
magnetic field B we find that the gap A x is vanishing as 



4A^x c (B) - x, 



(26) 



where x c (B) is defined by Eq. ( |21) , [The same result also 
holds for the energy minimum A y of the collective mode 
described by Imx yy (k, to; B).] Fig. [5] depicts the gap A x , 
Eq. ( 26 1, as a function of magnetic field for different coupling 
ratios x in the SOS phase. 



VI. DISCUSSION 



A. Summary 



In this work we have refined the theory, developed in 
Ref. Q3D of the QCP in a spin-orbital Hamiltonian for the A- 
site spinel compound FeSc2S4. The model exhibits an inter- 
esting quantum critical point: on increasing the second neigh- 
bor spin-exchange interaction J 2 it passes through a zero tem- 
perature phase transition from a spin-orbital singlet state to a 
magnetically and orbitally ordered phase. First, we consid- 
ered the effects of a weak nearest-neighbor exchange inter- 
action J 1, which induced a narrow region of incommensurate 
phase near the QCP. We studied the associated commensurate- 
incommensurate transition. Next, we included the effects of 
an external magnetic field. While the quantum critical point 
studied here seems similar to the one found in, e.g., bilayer 
Heisenberg antiferromagnets, its behavior under an external 
magnetic field is quite different. Namely, we found that a 
magnetic field suppresses magnetic and orbital order, and a 
transition from the ordered state to the spin-orbital singlet 
phase occurs at some critical field strength (see Fig. [2j. From 
these findings, we conclude that FeSc2S4, which is close to 
the quantum critical point, but in the spin-orbital singlet phase, 
does not show any field induced transition to an ordered state. 
Indeed, recent NMR experiments in fields up to 8.5T showed 
no signs of magnetic ordering. 10 Furthermore, we computed 
the dynamical spin susceptibility in the SOS phase by means 
of a random phase approximation. Averaging our results over 
the angular components of the wavevector we performed a 
comparison with available neutron scattering data on poly- 
crystalline FeSc 2 S4 samples and found reasonable agreement 
(see Fig. [4}. 
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B. Experiments 

1. Magnetic probes 

The theory expoused in this paper and Ref. 19 is broadly 
consistent with the results of a variety of magnetic probes 
applied to FeSc 2 S4. It explains the small but non-zero spin 
gap measured in inelastic neutron scattering and NMR 1/Ti 
relaxation rate measurements, as well as the temperature de- 
pendence of the uniform magnetic susceptibility. The present 
calculation of the dynamical spin susceptibility matches rea- 
sonably well with experiment. The lack of field-induced mag- 
netic ordering is also in agreement with the calculations in this 
paper. 

2. Specific heat and disorder 

The specific heat data on FeSc2S4 reveal several energy 
scales. The magnetic specific heat divided by temperature, 
C m /T, exhibits a peak at T » 6K. The integral of C m /T 
exceeds the spin-only entropy i?ln5, approaching instead 
i?(ln5 + In 2) for T > 60K, evidencing the 2-fold or- 
bital contribution. This is quite consistent with the present 
model. However, the lower temperature behavior is more 
complex. For T < 2K, experiments are fitted approximately 
by C m ~ AT + BT 2 - 5 , with the linear term dominant for 
T < 0.2K. The latter behavior appears at odds with the indi- 
cations of an energy gap of 1 — 2K in neutron scattering and 
NMR experiments. 

A possibility reconciliation of these observations is in the 
effects of disorder. Microscopically, we expect the dominant 
type of disorder to be inversion defects, in which the A and 
B sublattice atoms are interchanged. Inversion is very com- 
mon in spinels. To understand the effects of such defects, we 
apply general arguments based on the Landau expansion and 
the theory of disordered systems. These arguments depend 
very little upon the specific nature of the defects, other than 
that they are random, not very correlated, and do not break 
time-reversal symmetry. 

These conditions lead to an important observation: since 
the order parameters 1/^,0 are odd under time-reversal, dis- 
order couples only quadratically to them. Thus impurities 
behave, from the point of view of critical behavior, as ran- 
dom bonds rather than random fields. In three dimensions, 
it is known that in this case both phases are perturbatively 
stable to weak impurities. However, sufficiently close to the 
QCP, even weak impurities become non-perturbative. More 
formally, random bond disorder is a relevant variable at the 
QCR Physically, the most important effect of disorder is to 
locally break the degeneracy of the different ground states of 
the clean system within the ordered phase. For instance, a spe- 
cific impurity configuration might favor the p = 27r(l, 0, 0) 
state in one region and p = 27r(0, 1,0) state in another. Far 
from the QCP, the surface energy cost to create a domain wall 
between the two states overwhelms the random energy gain, 
and the system remains uniform. However, close to the QCP, 
the surface tension becomes small, and one expects the sys- 



tem to break into domains. Thus impurities induce a non- 
uniform disordered magnetic state, a "cluster spin glass", near 
the QCP. We expect, moreover, that this cluster spin glass state 
extends slightly past the QCP into the region of the SOS state 
of the clean system. This occurs because the system lowers its 
energy slightly more than in the clean case, by taking advan- 
tage of the impurities locally. 

This scenario provides a possible explanation of the spe- 
cific heat data. At low temperature, a T-linear specific heat is 
a generic feature of spin glasses. It should occur with a small 
coefficient A when disorder is weak. At higher temperature, 
one recovers approximately the intrinsic bulk clean behavior, 
which would be of the form C m — BT 3 f(A/k B T), where 
A is the energy gap, and f(5) is a monotonic scaling function 
satisfying f(0) = 1 and f(S) ~ S 7 / 2 e- s for S > 1. The T 3 
dependence is characteristic of the linearly-dispersion modes 
at the QCP, which is cut off by the gap. It seems plausible 
that the experimental observed T 2 5 dependence reflects the 
attempt to fit such a form to a single power-law. If the im- 
purities are not too weak, it is also plausible that they modify 
the T 3 behavior somewhat. In any case, the overall behavior 
seems reasonably in line with theoretical expectations. 



C. Directions for future work 

The theory in this paper (and Ref. [191 appears to give a con- 
sistent explanation for the experimental results on FeSc2S4. 
However, there are a number of directions that could be ex- 
plored in the future. It would be desirable to have a di- 
rect proof of the postulated spin-orbital entanglement in the 
ground state of FeSc2S4. Theoretical proposals and experi- 
mental studies to this end would be welcome. Given the small- 
ness of the gap in FeSc2S4, there is a possibility that it might 
be driven across the QCP by pressure, which would be very 
exciting. Looking more broadly, it appears that the mecha- 
nism for quantum criticality described here could apply at the 
very least to any material with Fe 2+ ions in a tetrahedral en- 
vironment. It would be interesting to survey such compounds 
for signs of this physics. 



Acknowledgments 

This work was supported by the DOE through Basic En- 
ergy Sciences grant DE-FG02-08ER46524. LB's research fa- 
cilities at the KITP were supported by the National Science 
Foundation grant NSF PHY-0551 164. 



APPENDIX A: HUBBARD-STRATONOVICH 
TRANSFORMATION AND LANDAU ACTION 

In this appendix, we derive the effective Landau action, 
Eqs. Q and |6]), from the microscopic Hamiltonian (JT| using 
the Hubbard-Stratonovich method to decouple the exchange 
interactions. In this way, we can relate the coefficients in 
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Eqs. Q and |6]) to the microscopic exchange coupling pa- 
rameters Ji and the spin-orbit coupling A. We consider the 
J1-J2-X model on the diamond lattice, 

W = 53JiSi-S 3 -+ 2 + (Al) 



(ij) 



am 



where the brackets (ij) and ((ij)) denote the summation over 
first and second nearest neighbors, respectively. The on-site 
spin-orbit coupling term Hq is given by Eq. ( fTcj i. The partition 
function reads 



Z = Tr exp 



(A2) 



respectively. We decouple the exchange interaction by intro- 
ducing the auxiliary field 4>i and transforming the partition 
function to 



Z = [v<fK$ J **-+'T*e-fiXiW>-*i-*i). (A3) 



Expanding around the saddle point yields 



£></>exp [-S e ff] , 



(A4) 



with the exchange coupling matrix Jy. Here, Jij — J\ or 

J2, when ij connects first neighbor or second neighbor sites, with the effective action 



Seff = I dr 
Jo 



2 J l3 <t>i-<t>i + x , A ' 2 



A 3 



Assuming J\ <C J2 and expressing the <fi^ field in terms of the staggered magnetizations ■0 Mj a 



we obtain the quadratic part of the free energy density 



(A5) 



(A6) 



/ (2) (M) -EE(^-5 + 3^(1 + §2) k l) [l^A,a(k)| 2 + |^B, Q (k)| 2 ] + ^fe ^,a(-k) • ^ B)0 (k). (A7) 



k a 



We note that Eq. \A1\ is compatible with Eqs. Q and ([6j, 
which we derived using symmetry considerations. By com- 
paring the coefficients in Eqs. Q and |6]) to those in Eq. \A.1\ 
we find the following relations 



r 
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1/(32J 2 )(1+£|), 
l/(8J 2 )-2/A, 

Ji/(16J|); 



(A8) 



Similarly, expressing the quartic terms in Eq. ( |A5| l in terms 
of the staggered magnetizations one can show that the coeffi- 
cients <?i and g-2 are given by 



.91 = 4/A 3 , 
92 = -2/A 3 . 



(A9) 
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